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We use the QCD sum rules to analize the hadronic decay 73^(2317) — ■> D^-k , in the hypotesis 
that the 1)^(2317) can be identified as a four-quark state. We use a diquak-antidiquark current 
and work to the order of m s in full QCD, without relying on l/m c expansion. We find that the 
partial decay width of the hadronic isospin violating mode is proportional to the isovector quark 
condensate, (0\dd — uu\0). The estimated partial decay width is of the order of 6 keV. 

PACS numbers: 11.55.Hx, 12.38. Lg , 13.25.-k 



The strange-charmed mesons L»+ 7 (2317) and L>+ 7 (2460) with spin parity 0+, 1+ [J H H are lighter than the 
prediction of the very successful quark model for the charmed mesons . One interpretation is that this is evidence 
for a chiral symmetry where the mass gap between the 0~ and 1~ states equates the mass gap between the + and 1 + 
states @ . Other interpretations are related with the underlyingstructure of these mesons, which has been extensively 
debated. They have been interpreted as conventional cs states rZk B , u i llCl fill fl2l Ha. ll4L Ha . Ilq . two-meson molecular 



O 



X 



state |17lll8| . D - K- mixing |19| or four-quark states 12^121112^12^12^,1251,12^,1271]. 

Because of their low masses, these two states are lower than the DK and D*K thresholds. Therefore, their strong 
. decays must proceed th roug h isospin violating effects. There have been some discussions of their decays within the 
quark model @, [H Q El H3 and QCD sum rules 01- I n & ll these studies but the isospin violating effects 
were considered through the rj — ir° mixing. However, if these mesons are considered as four-quark states, in a QCD 
sum rule calculation, the isospin violating effects can be introduced through the mass and quark condensate difference 
, between the u and d quarks. 

In a recent calculation [2j| the scalar- isoscalar meson D^ / (2317) were considered as a S- wave bound state of a 
diquark-antidiquark pair. As suggested in ref. [2S|. the diquark was taken to be a spin zero colour anti-triplet. The 
corresponding interpolating field is: 



js = eabc ^ ec [{ulC^c h ){u d ^Csl) + {dlC^c b ){d d ^Csl )] , (1) 

o 1 

where a, b, c, ... are colour indices and C is the charge conjugation matrix. In ref. 23], using the QCD sum rule 
(QCDSR) formalism [2^. l30 H3l| . it was shown that it is possible to reproduce the experimental mass of the meson 
Z?j"j(2317) using this four-quark state picture. Here, we extend the calculation done in ref. 0] to study the vertex 
associated with the decay Z?^ 7 (2317) — > Dfir . 

The QCDSR calculation for the vertex, D^j (2317) ir° , centers around the three-point function given by 



r„(p,p',q) = [ d*xd A y e*'* (0\T[j Ds ( x )jf,(y)jl (0)]|0>, (2) 



where p = p' + q and the interpolating fields for the pion and D s mesons are given by: 
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■?5m = —^( U altil5U a ~ <4 7^75^4), JD, = is a "/5C a . (3) 



The fundamental assumption of the QCD sum rule approach is the principle of duality. Specifically, we assume 
that there is an interval over which the above vertex function may be equivalently described at both, the quark level 
and at the hadron level. Therefore, the underlying procedure of the QCD sum rule technique is the following: on 
one hand we calculate the vertex function at the quark level in terms of quark and gluon fields. On the other hand, 
the vertex function is calculated at the hadronic level introducing hadron characteristics such as masses and coupling 
constants. At the quark level the complex structure of the QCD vacuum leads us to employ the Wilson's operator 
product expansion (OPE) . The calculation of the phenomenological side proceeds by inserting intermediate states for 
D s , 7T° and D s j, and by using the definitions: 

2 p 

<0|J5>°(9)> = ig^, (0\ jDs \D s (p>)) = 2°+^ (°\jD SJ \D sJ (p)) = As. (4) 
We obtain the following relation: 

r^ hen (p,p' , q) — — Ds ^ Da — - 9 Ds jDa7r _j_ continuum contribution , (5) 

(m c + m s )(p 2 - m 2 Dsj )(p' - m 2 Dg )(q 2 - m%) 
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where the coupling constant gD SJ D s -w is defined by the on-mass-shell matrix element 

(D s ir\D sJ ) = Qd s jD 3 7r ■ 



(0) 



The continuum contribution in Eq.Q contains the contributions of all possible excited states. 

For the light scalar mesons, considered as diquark-antidiquark states, the study of their vertices functions using the 
QCD sum rule approach at the pion pole |T TL I30L 133, l33| . was done in ref.[34|. In Table I we show the results obtained 
for the different vertices studied in ref. |34| . as well as the experimental values. 

Table I: Numerical results for the coupling constants 



vertex 


g(GeV) 


g exp (GeV) 


<jir + n~ 


3.1 ±0.5 


2.6 ±0.2 




3.6 ±0.3 


4.5 ±0.4 


foK+K- 


1.6 ±0.1 




/ 7r + 7r~ 


0.47 ±0.05 


1.6 ±0.8 



From Table I we see that, although not exactly in between the experimental error bars, the hadronic couplings 
determined from the QCD sum rule calculation are consistent with existing experimental data. The biggest discrepancy 
is for g/oTr+Tr- and this can be understood since the /o — > tt + tt~ decay is mediated by gluon exchange and, therefore, 
probably in this case a, s corrections, which were not considered, could play an important role. 

Here, we follow ref. |34| and work at the pion pole. The main reason for working at the pion pole is that the matrix 
element in Eq.@ defines the coupling constant only at the pion pole. For q 2 ^ one would have to replace gD s .,D 3 n, 
in Eq.©, by the form factor gD SJ D s Tr(g 2 ) and, therefore, one would have to deal with the complications associated 
with the extrapolation of the form factor 35, 36]. The pion pole method consists in neglecting the pion mass in the 
denominator of Eq. (JSJ and working at q 2 = 0. In the OPE side one singles out the leading terms in the operator 
product expansion of Eq.@ that match the l/q 2 term. In the phenomenological side, in the structure |f we get: 



T phen (p 2 ,p' 2 ) = 



s)(p 2 



'D 



+ 



Pcont(p 2 , U) 



du. 



u — p 



(7) 



In Eq.lJTJ, p C ont{p 2 ', u), gives the continuum contributions, which can be parametrized as p C ont{p 2 , u) = J^pH 0(u — uq) 
|37ll38| . with s an d uq being the continuum thresholds for D s j and D s respectively. Since we are working at q 2 = 0, 



we take the limit p 2 = p' and we apply the Borel transformation to p 2 — > M 2 to get: 



YP hen (M 2 



\ s m 2 D 



fD s Fn gD s 



(m c + m s )(m 2 D 
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>j' m ) + A e~ s ^ M2 + / p cc (u) e- ulM *du 



where 



(s Q - u) 



du, and p cc (u) 



(so - u) ' 



(8) 



(9) 



stands for the pole-continuum transitions and pure continuum contributions. For simplicity, one assumes that the pure 
continuum contribution to the spectral density, p cc (u), is given by the result obtained in the OPE side. Asymptotic 
freedom ensures that equivalence for sufficiently large u. Therefore, one uses the ansatz: p cc (u) — Pope{ u )- In Eq.JSJ, 
A is a parameter which, together with gD SJ D a ir, have to be determined by the sum rule. 

In the OPE side we work at leading order and deal with the strange quark as a light one and consider the diagrams 
up to order m s . To keep the charm quark mass finite, we use the momentum-space expression for the charm quark 
propagator. We calculate the light quark part of the correlation function in the coordinate-space, which is then 
Fourier transformed to the momentum space in D dimensions. The resulting light-quark part is combined with the 
charm-quark part before it is dimensionally regularized at D = 4. Singling out the leading terms proportional to 
Qii/g 2 : we can write the Borel transform of the correlation function in the OPE side in terms of a dispersion relation: 



r 



OPE 



(M 2 ) = 



Pope 



(tt) e' u/M2 du , 



(10) 
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where the spectral density, pope, is given by the imaginary part of the correlation function. Transferring the pure 
continuum contribution to the OPE side, the sum rule for the coupling constant, up to dimension 7, is given by: 



1 



2 4 ^ 2 



du e- u > M2 u 1 



m c m s 
' 8vr 2 



du e -u/M> h _ m* 



m c Q c -,n 2 JM 2 



M 2 



with 



C = 



9D s .,D s tt- 



(m c + m s )(m 2 Dsj - mjj j 

In Eq. (|ll() . 7 measures the isospin symmetry breaking in the quark condensate: 

_ (0\dd-uu\0) 
1 ~ (0|tm|0) ' 



(11) 



(12) 



(13) 



and is the source of nonperturbative isospin violation in the OPE of the vertex function. The value of 7 has been 
estimated in a variety of approaches [32U391 40 1 . with results varying over almost one order of magnitude: — 1 x 10~ 2 < 
7 < —2 x 10~ 3 . A more recent calculation [4l| gives a bigger (in module) value: 7 = —2.6 x 10 -2 . 

In the numerical analysis of the sum rules, the values used for the meson masses, quark masses and condensates 



are: m DsJ = 2.317 GeV, m D , = 1-968 GeV, m s = 0.13 GeV, m c = 1.2 GeV, (qq) 



-(0.23) 3 GeV , (ss) = 0.8(qq). 



For the meson decay constants we use — y/2 93MeV and fo 3 = 0.22 GeV (obtained using uq = 6 GeV 2 |42|1. For 
the current meson coupling, Ag, defined in Eq.Q we are going to use the result obtained from the two-point function 
in ref. 23]. Considering 2.6 < ^ < 2.8 GeV we get A s = (2.9 ± 0.3) x 10" 3 GeV 5 . 

In Fig. 1 we show, through the dots, the right-hand side (RHS) of Ea. fTT|l . for u = 6 GeV 2 and 7 = — 1 x 10~ 2 
40], as a function of the Borel mass. We use the same Borcl window as defined in ref. 231. 
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FIG. 1: Dots: the RHS of Eg. illlll . as a function of the Borel mass. The solid line gives the fit of the QCDSR results through 
the LHS of Eq.{nj- 

To determine Qd s jD s 7t we fit the QCDSR results with the analytical expression in the left-hand side (LHS) of 
Eq.JITJl: 

C (e- m °J M2 - e- m °sJ /M *) + A e-WM 2 ( (14) 

and we get (using ^ = 2.7 GeV): C = 8.34 x 10~ 6 GeV 7 and A = 6.96 x 10~ 6 GeV 7 . Using the definition of C 
in Eq.(H2I) and X s = 2.9 x 10 -3 GeV 5 (the value obtained for = 2.7 GeV) we get 5d 3J d s tt = 50 MeV. Allowing 
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so to vary in the interval 2.6 < y/s~o < 2.8 GeV, the corresponding varyation obtained for the coupling constant is 
45 MeV < g DsJ D s TT < 55 MeV. 

Fixing y/sTj = 2.7 GeV and varying the quark condensate, the charm quark and the strange quark masses in the 
intervals: -(0.24) 3 < (qq) < -(0.22) 3 GeV 3 , 1.1 < m c < 1.3 GeV and 0.11 < m s < 0.15 GeV, we get results for 
the coupling constant still between the lower and upper limits given above. However, varying the value of 7 form 
— 1 x 10~ 2 to the more recent value given in |4l|: 7 = —2.6 x 10~ 2 , and keeping the other parameters fixed we get 
9d,jD s tt = 130 MeV. On the other hand, if we use the smallest (in module) value allowed for 7: 7 = —2. x 10~ 3 , we 
get gD SJ D s -K — 10 MeV. Therefore, the biggest source of uncertainty in our calculation is the value of 7. In all cases 
considered here, the quality of the fit between the LHS and the RHS of Eq.^TJ is similar to the one shown in Fig.l. 

The coupling constant, gD SJ D s m is related with the partial decay width through the relation: 

r(D+ (2317) -> Dfn°) = \ g 2 Ds J\(m 2 Dsj , m 2 Ds , m 2 ), (15) 

where X(a, b, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 2bc. Considering 7 = — 1 x 10~ 2 , which was the value found in j4jj to 
be consistent with the neutron-proton mass difference in a QCDSR calculation, and allowing (qq), m c , m s and sq to 
vary in the ranges discussed above we get: 

r(£>+ 7 (2317) -> DfiT ) = (6±2)keV. (16) 

However, it is important to state that, if the value for 7 found in ref.|4l| proves to be correct, then the partial decay 
width could be as large as r(Z)jj(2317) — * Dfn°) ~ 40keV, in agreement with the QCDSR calculation done in 
ref.^(|, where the meson Dfj (2317) is considered as a ordinary cs state. 

In Table II we show the partial decay width obtained by different theoretical groups. The first five calculations 
assume a cs picture for Dfj (2317), while the last two assume a four-quark picture for it. 

Table II: The decay width r(D+ J (2317) -> D+n°) from various theoretical approaches. 



ref. 


m 


\m 


[14] 


[15] 


\M 


W 


this work 


T (keV) 


21.5 


7± 1 


- 10 


16 


39 ±5 


10 - 100 


6±2 



From the results in Table II we see that we can not get a definitive answer about the structure of the D s j (2317) 
meson from its strong decay width, since in both pictures: ordinary cs or four-quark states, different approachs can 
give results varying from a few keV to a hundred keV. 

We have presented a QCD sum rule study of the vertex function associated with the strong decay Djj(2317) — > 
D+n°, where the D+ 7 (2317) meson was considered as diquark-antidiquark state. We found that the source of isospin 
violation in our calculation is the parameter 7 = ^ofsJ l oj ^ ' which measures the isospin symmetry breaking in the 
quark condensate. Since, in our approach, the partial decay width is directly proportional to 7 2 , and since there is a 
large uncertainty in the value of 7, considering 7 in the range —2.6 x 10 -2 < 7 < —2 x 10~ 3 we get the partial decay 
width in the range 0.2 keV < r(D+ J (2317) -> Dfir ) < 40keV. However, from other QCDSR calculation, we believe 
that the value of 7 is ~ —1 x 10~ 2 , which gives the result shown in Eq. l|16|l . 

As a final remark we would like to point out that if, instead of using a isoscalar current, we have used a isovector 
current for D s j(2317) (as suggested in ref. the difference in Eq.(^TJ would be a factor 2 in the place of 7. In 

this case the decay would be isospin allowed and the partial width would be ~ 230 MeV, much bigger than the 
experimental upper limit to the total width r ~ 5 MeV. 
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